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ABSTRACT 


A combined analytical-numerical procedure for determining 
precise mean orbital elements is presented and applied to the 
orbits of GEOS-I and GEOS-II. The precision of the mean 
semi-major axes of these orbits is a few tens of centimeters 
when optical flash data is used to determine 2 day orbital 
arcs. Four day Minitrack orbits give mean semi-major axes 
of a few meters precision. The mean orientation parameters 
(i, ft) are obtained to a precision of about O'.'l (~3m) or 
better from the optical orbits. This precision is adequate 
for determinations of tidal parameters, particularly in the 
case of GEOS-II where the tidal perturbation of the inclin- 
ation is 10". 
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1 . Introduction 


The most familiar method for study of the long periodic 
and secular perturbations of an orbit is to compare the 
variation of the mean elements of the orbit with the changes 
predicted by theory . This method has been used with success 
by many investigators to determine geopotential coefficients, 
and atmospheric and tidal parameters. It is the purpose of 
this paper to consider the problem of transforming osculating 
elements to mean elements with minimum loss of accuracy. 

We begin by considering how accurately the osculating 
elements of an orbit can be determined. This is a complex 
question depending on the distribution and number of data, 
orbital arc length, model parameters used, etc., but some 
generalizations are po-ssible due to recent efforts at 
Goddard Space Flight Center. 

A recent paper by Marsh and Douglas C 1 9 7 1 ) concerning 
GEOS- I and GEOS- II shows that for orbital arcs of a few days 
duration, uncertainty of the geopotential, particularly 
resonant coefficients, is the most important error source. 

If resonant coefficients are adjusted the orbital error 
along-track can be reduced to a few milliseconds in time for 
a 5-1/2 day orbital arc heavily observed by optical trackers. 
This is equivalent to about 15 meters. Gaposchkin and Lambeck 
(1970) report position accuracy of similar magnitude for 
near-Earth orbits. 

Marsh and Douglas (1971) show further that the orbit error 
tends to be of high frequency. The dominant period of the 
error is the period of the orbit itself. Thus a least squares 
fitting procedure will tend to treat this error as "noise" 
in the sense that the satellite is ahead in its orbit as 
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often as it is behind its "true" position. We should expect the 
orbital elements to be very precise (if not accurate), 
indeed much more so than the satellite position itself. 

One should anticipate obtaining mean elements from relatively 
short orbital arcs that are precise to a few meters or less. 
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2. The Definition of the Mean Elements of an Orbit 

We wish to remove the high frequency perturbations from 
a sequence of osculating element sets in order to study the 
long period and secular variations of the orbit. This 
will involve removal of the short periodic effects, i.e., 
those with periods equal to or less than the orbital period, 
and effects introduced by the rotation of the Earth, the 
m-daily effects of tesseral harmonics. 

The dominant short periodic effects are due to the second 
zonal harmonic. The amplitude of these effects on satellites 
such as GEOS- 1 and GEOS- 1 1 is many kilometers. The first 
order effects are easy to remove by using the equations of 
Brouwer (1959), or, more generally, by using those of 
Kaula (1966) . The first order short periodic and m-daily 
effects of tesseral harmonics are also easy to remove 
analytically, particularly if the previously mentioned 
development of Kaula (1966) is used. However, short periodic 
effects of the sun and moon, drag, radiation pressure, second 
order effects of oblateness, and the interaction of oblate- 
ness with other perturbations must also be considered. 
Obviously, the analytic calculation of all of these small 
effects is complex, particularly if accuracy at the 1 meter 
level or better is. required. Thus we chose to use analytic 
techniques to remove only the dominant oblateness and 
tesseral harmonic perturbations, and to employ a numerical 
method for removal of the other perturbations. Although a 
purely numerical method to remove all high frequency per- 
turbations may be theoretically possible, we shall see below 
that very great efficiency and accuracy is obtained by the 
combined method. 

It is common to think of the mean elements of an orbit 
as "averages" in some sense. However, examination of the 
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gravitational disturbing function shows that the perturbations 
do not all average out over the same period. For example, 
during an orbital revolution, the sun and moon move and the 
Earth rotates significantly. The time over which the short- 
periodic perturbations average out is slightly different 
than a revolution, and most importantly, is different for 
each perturbing source. Moreover, in satellite theories 
we usually take into account the motion of the node and perigee 
in the computation of short-periodic terms, that is,, 
solutions take the form of forced oscillations about a 
secularly precessing Kepler ellipse. Thus the high frequency 

perturbations of the geopotential have the frequencies 
(Kaula , 1966) : 


iu + jM + k(ft - 0) 2.1' 

where i, j and k are integers, 0 is the rotation rate of the 
Earth and <L, M, fi are the Kepler element rates. Neither 
the short-periodic and m-daily geopotential perturbations 
average out over any common period. In the language of 
electrical engineering, we really need to filter the oscu- 
lating elements with an ideal low-pass filter. Removal of 
high frequency terms by very accurate analytic methods approx- 
imates such a filter. A purely numerical averaging filter 
has relatively poor characteristics because of the lack of 
any unique period over which all frequencies exactly average. 
However, by confining the numerical averaging to small (<50m) 
effects, the error introduced by the averaging is tolerably 
small . 

Considering these remarks, our scheme for determining mean 
elements takes the following form: 

We first generate an ephemeris in terms of osculating 
elements at one minute intervals for 1 day. From each set of 
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these elements are then subtracted the short periodic 

oblateness, m-daily, and resonant perturbations. These 

preliminary mean elements show a variation of 30-50 meters 

for the GEOS satellites. The preliminary elements are then 

fitted by least squares to a secularly precessing Kepler 

_ _ _ • • 

ellipse so that the nine parameters a,e, I, w 0 , to, ft, ft, 

M q , and M are determined. The epoch of these elements is 
taken to be the mid-point of the averaging interval; 
of course the rates w, M and ft are used to transform w, 

M and ft to this time. In this way long periodic and secular 
variations are properly represented in the averaged 
elements. 

The necessity for this combined scheme is shown in 
Figure 1. Mean semi-major axes calculated by purely numer- 
ical averaging (X) and the combined method (•) are shown for 
2 day GEOS- II optical data arcs. In the elements obtained 
purely numerically the subtle decay of the semi-major axis 
is not even detectible. 

Figure 2 shows the mean semi-major axes of GEOS- I during 
1965-66. Note that the precision is about 25 cm. GEOS-I 
suffers very large radiation pressure perturbations, as- is 
obvious from the increase of more than 20m in the semi- 
major axis in early 1966. Figure 2 has been very useful for 
geodetic investigations involving GEOS-I because inconsistent 
data arcs are clearly distinguishable (for example, refer to 
the outlying arcs in April 1966). 

Figure 3 shows the mean semi-major axes of GEOS-II for 
1968 determined from 2 day optical data arcs. All arcs in 
this paper were reduced using BIH Polar Motion and UT1 
data, the 1969 SAO Standard Earth (Gaposchkin and Lambeck, 
1970) gravity model, and a worldwide network of SAO, NASA, 
and International participants optical tracking stations at 


5 



coordinates estimated by Marsh, Douglas and Klosko (1971). 

Note that the orbit of GEOS- 1 1 is much more stable than that 
of GEOS- I against radiation pressure perturbations, although 
the decay in semi-major axis is highly variable. 

Figure 4 presents the mean semi-major axes of GEOS-II 
where optical data was available in 1969. The precision is 
poorer than in 1968. No explanation is available. 

It would be tempting to conclude from Figures 1-4 
that the resonant coefficients for GEOS-I (12th order) and 
GEOS-II (13th order) must be known to extreme accuracy , 
because their effects on the mean elements were essentially 
perfectly removed. However, the beat periods for these orbits 
are only about 7 days, and the effect of an error in 
the coefficients tends to be reduced because of the relatively 
long averaging time. The recent investigation by Marsh 
and Douglas (1971) shows that the SAO 1969 Standard Earth 
models about 90% of the resonance effect for GEOS-II. 

The remaining uncertainty should be detectible in the GEOS-II 
mean semi-major axes, but the smoothing, procedure has 
obliterated the effect. 

Figure 5 shows the semi-major axes for GEOS-II in 1970 
obtained from 4 day Minitrack-determined orbits. The 
scatter is about 2m, a precision sufficiently accurate for 
studies of atmospheric density. The orientation elements are 
less well-determined from the Minitrack data (i.e. about 10 
arc seconds) . 

Mean elements for GEOS-I and GEOS-II obtained from 
optical arcs are given in Tables 1-4. The Minitrack 
elements are given in Table 5. All mean elements are referred 
to the true equator and equinox of date. 
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Figure 1 . GEOS-II Mean Semi-Major Axis 
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MOOIFIED JULIAN DAY 
EARTH RADII (a„ = 6378155 ml. 


Table 2. GEOS-II Mean Elements from 2-Day Optical Data Arcs in 1968 
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A 

E 

INCL 

0*EG A 

NODE 

PEAK 
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1 05*709563 
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123.119100 
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357. 1 50473 
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105.732057 

24.276055 
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112.016673 


1 .206061 99 
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331 .419661 
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•C3181 2143 
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1 • 206061 69 
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26.157067 
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i .208061 er 
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333.941722 
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Table 3 . GEOS-II Mean Elements from 2 -Day Optical Data Arcs in 1969 
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Table 4. GEOS-II Mean Elements from 4-Day Minitrack Data Arcs in 1970 
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